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ABSTRACT 

As  an  extension  of  an  earlier  report  the  Green's 
function  Is  evaluated  for  the  Lundqulst  equations  linearized 
about  uniform  magnetic  field,  constant  matter  density  and 
zero  flow  velocity.   It  is  assumed  that  all  quantities  are 
functions  of  two  space  variables  and  time  only.   In  the 
general  magnetic  field  configuration  considered  here  a  pure 
Alfven  disturbance  no  longer  exists;  there  is  Instead  a  wave 
with  properties  of  both  the  Alfven  and  fast-slow  disturbance, 
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On  the  Green's  Function  for  Two-Dimenslonal 
Ma gnetohydro dynamic  Waves.  II 


In  a  previous  paper  v/e  studied  the  fundamental  solution  of 
the  linearized  Lundquist  equations  for  a  restricted  two  dimen- 
sional situation.  One  of  the  essential  simplifications  resulting 
from  the  assumptions  made  was  the  decoupling  of  the  Alfven  and 
fast-slow  (or  "longitudinal")  waves.  Here  we  wish  to  consider 
the  most  general  two  dimensional  situation  in  which  the  Alfven 
v^aves  are  coupled  to  the  fast-slow  waves.   For  the  "physical" 

part  of  the  fundamental  solution,  we  shall  again  give  an  algebraic 

2 
function  plus  standard  distributions.   We  shall  draw  heavily  on 

the  material  presented  in  I  and  omit  all  arguments  which  are 
parallels  or  direct  extensions  of  those  found  there. 

We  carry  over  all  the  hypotheses  of  I  but  one,  we  now  assume 
that  the  constant  magnetic  field  need  not  lie  in  the  x  y  plane. 
Instead  (and  with  somewhat  irregular  notation)  we  assume  (compare 
1(2)) 

(1)       B^  =  B^(1,0,€)  =  B^. 

With  this  extension  we  now  have  the  general  two-dimensional  situa- 
tion. There  is  one  other  special  two  dimensional  situation,  where 
B  =  (0,0,B(x,y)),  which  has  been  studied  even  for  the  non-linear 
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equations-^  and  which  we  ignore.  Besides  the  symbols  in  I  we 
need  one  new  speed  similar  to  the  Alfven  speed: 


2   ^o  ^^  +  ^^ 
(2)       a2=^ 


The  equations  of  motion  are  just  the  same  as  1(1)  and  the 
equations  for  the  transforms  are  the  same  as  l(5)i  they  are: 

(3a)  cop   -    (k'u)   =  iSp 

(3t))  CDU  -  a^  Icp  +  A^  B'x(Bxk)   =  iS^ 

(5c)  a^  +  Tc  x(uxb)   =  iSg   . 

With  b  given  its  new  meaning  by  (l),  we  may  no  longer  split 
off  the  z  components  of  u  and  B  as  in  I.   Instead  we 
use  a  decomposition  procedure  essentially  equivalent  to  that 
described  by  H.  Grad. 

From  (3)  we  derive 

(4a)      co(k.u)  -  al   k^p  +  A^(  (B.k)(b.k)-  k^(^.B))  =  l!?.S^ 
(4b)      ay{t'^   -   a^(k.'b)p  =  iB'.S^ 

(4c)      od(b."^)  +  (u."5^(k.'b)  -  b^(k.^  =  iB'.Sg 

(4d)      cja(I^,B)  =  i^-^, 

and  with  (3a)  we  have  a  closed  system  for  p,    (^-u),  and  (^-B), 
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u 


so  that 


(5a)    p  =  i  (o^Cco^-  Afk^)Sp+  co2(k.Sy)-  k2A^(b.k)(b.S^) 


,2, 


+  aiAj(kKb.kxSg)) 

(5b)    (u.b)  =  i  (a)(co2-  Ajk2)(b.S^)+  a^a)(bxk.kxs^) 

+  k.b(a3^-  k^A^)S  +^^A^  b.k(k  x^-kx  ^) ) 

(5c)    (B'b)  =  4  (03(0)2-  k2a^)tp.b  +  o^S  al{k%^- .{k-h)^) 

+  0)  (bxk«bxs^)+  ajb«k(k^b«kx  s  )) 


+  ^  (k.S3)(^.!f)A2  (a^C^.k)^-  bV) 


where 


(6)     D  =  D(o),k)  =  cJ^-   o.2(a2  +  A2)k2  +  l^^i^^^fa^Al 


(Compare  with  D  as  defined  In  I. )  In  (5c)  the  second  term 
Is  the  "unphyslcal"  term.  We  may  evaluate  all  the  transforms 
in  (5)  directly  by  the  methods  of  I. 

If  we  can  give  uxb  and  Bx.b,  possibly  in  terms  of 
p,  (u*^),  {B'h) ,  then  we  do  have  p,  u  and  B,  but  again 
from  (3)  we  derive 
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(7a)    ai(uxb)  +  A2(Bxb)(b.k)  =  i{S^Kb)+  k><b(a^p  +  A^(b'B)) 


(7b)    ui{Bxt)   +   (uxb)(b-k)  =  i(S3Xb). 

It  Is  possible  to  give  many  equivalent  representations  of 
the  solutions  of  (7).   We  choose  an  unsymmetrlcal  one  which 
emphasizes  certain  properties  of  the  solution.   In  the  form 
chosen  we  see  easily  that  for  e  =  0  the  results  coincide  with 
I  and  thus  we  have  separated  the  part  of  the  Alfven  wave  that 
does  not  vanish  as  e  — >  0.  We  define  the  vectors 

(8a)  t  =   (0,1,0) 

(8b)  n  =  (-€,0,1), 


-A  _x  -A  _X  _i 

so  that  b,J,n  form  a  right  handed  orthogonal  set.   With   (u'b) 

and   (B'b)  given  we  need  only  specify  the  quantities   (u-j), 

(u-n),  {B't),    (B-n).  With  the  definition 


(9)  A  =  0.2-  aI    (i?.^)2 

we  find 

(10a)   (^-t)  =  4  |a^co((k.S^)  +  (DSp)(k.j')  +  0.(0)^-  aV)(S^- j") 

+  A^  a^(k.b)(k.n)(b.kxag)+(n.kxag)A^(oD^-a^(k.b)2) 

+  —   (n.k)A^cc^  ja3(b.kXS^)-A^(if-b)(kxSg.b)  V 
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(10b)      (^.n)    =  i[(t-l^x^g)A2  +  S^.^j  +  i|  f{h^-l){t^^t^)a.^ 
-  aQ(j-k)(n"k)(n"kxS3)  +  a^a)(k«n)  (k-S^  +  a)Sp)f 

+  ^^  I  (b2-l)(i^.^)2(^.i^x:^)-b2(^.lf)(t.k)(n.I^>.  Sg) 


+  CD(k-n)(b'>ck''b  X  s'  ) 


(10c)      (B'J^)   =  I     -a^(5^"i^)(l^.j^)(i?.^^4<oSp)4<D(^.3')(a)2-aV) 


-(S^-f)(b-k)(co2.aV) 


la3^A^(k»n) 


AD 


tD(Sg>^k-'5')H-(b"k)(b-kKSy) 


+  ^  (k-tB)|A^(k.t)(a2(I^.^)2-  co^b^) 
(lOd)      (B-n)   =  M<D(S3-n)-(b-k)(S^-n)j  -i  (k-n)  (k-b)a^(k-S^+a)Sp) 


p      p 

+  ^!L^  j  »(k.n)(k.b)(b<k.bxt^)4<o(b2-i)k2(S3-n) 


AD 


_\,  .A  _^,  ,_i  ->x 


+  a.(S„'J)(k-n)(k-j)'. + 


"B 


i(k-S^)(k-n)A 


3 


cdD 


°/a2(1^.^2_^2^i 


In  the   limit      e  — >  0,     k'n  ^   0     and     b   -   1  — >  0,      so  that 
the  terms  with  denominator     AD     cSisappear  and  the   solutions 
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go  over  to  those  In  I,  that  is  1(4)  and  1(6).   In  (10a)  and 
(10b)  only  the  divergence-free  part  of  Sg,  l^^  Sg,  appears. 
In  (10c) and  (lOd)we  have  split  the  expressions  into  physical 
and  "unphysical"  terms.   The  unphysical  ones  have  the  factor 
k'Sp.   In  (10c) the  remaining  physical  terms  involve  only  the 
combinations  ^^S^     or  ^'J   and  in  (lOd)on  rewriting 
{'^^'t){^'n)   =  (^gXl^-"5')  +  (lt.t)(^B-n)   only  T^x^^^  or  n'tg 
appear.  Thus  we  may  parallel  the  argviments  of  I  and  assert 
we  have  made  the  separation  into  physical  and  unphysical 
terms. 

We  may  evaluate  the  terms  in  the  transforms  with  the 
denominator  A  directly.   We  express  those  with  denominator 
D*  quite  simply  in  terms  of  the  functions  in  I.   If  we  set 


(lit)     tl   =  ili  -  I  /(a^  .  4f-  ,si  ,i 


then  D(cD,l<)  =  D{(a,'^,a    ,A  ),  where  D  refers  to  the  denomina- 
tor considered  in  I.   It  is  again  possible  to  show  that  the 
discontinuity  planes  induced  by  the  double  points  of  the 
normal  speed  locus  do  not  appear  and  we  may  carry  over  formulas 

1(11)  through  1(22)  with  the  replacement  (11),  and  thus  func- 

^    ^    ^ 

tions  G^^,  G^„,  G„„  are  defined.   The  singularity  surfaces 
J^-^        Ay    yy 

and  lacunas  for  this  case  given  by  I(19),  1(20)  after  the 
modification  (11);  thus  they  are  similar  to  those  in  I.  The 
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"unphyslcal"  terms  are  again  similar  to  those  In  Ij  as  before 
we  shall  not  Include  them  in  the  final  answers.  We  must  next 
consider  the  more  complicated  transforms  with  denominator  AD. 

The  methods  used  in  I  for  the  explicit  evaluation  of  a 
class  of  Fourier  transforms  required  that  the  functions  (or 
transforms)  depend  non-trivially  on  three  independent  variables. 
Although  the  Alfven  wave  depends  on  two  variables,  in  I  it 
was  independent  of  the  longitudinal  wave  so  we  could  determine 
it  separately.   Now  there  are  terms  in  the  Green's  function,  i.e., 
those  with  denominator  AD,  which  are  coupled  Alfven-longitudinal 
waves  and  we  need  some  method  of  handling  them.   We  shall  go 
through  a  circuitous  procedure  to  show  that  the  finite  part 
interpretation  of  a  singular  integral  is  again  appropriate. 
We  shall  proceed  to  make  a  series  of  formal  manipulations  of 
the  transforms  until  we  can  recognize  the  answer,  and  finally 
we  shall  recompose  the  answer. 

By  the  simple  identity 

P   ?  2    2   2 

Aq  ^        (-1)      ,    ^  -  ^o  ^ 

AD     (kxb-k><b)A   (kKb-kxb)D 

we  are  able  to  make  a  formal  decomposition  of  the  Alfven  and 
longitudinal  waves.   We  study  a  transform  of  the  form 

A^  0.2-^^  (k,)-(k  )P     j  ^        "^ 

(12a)   T  =  -2 —^ ^—  <  P<1 

AD 

L         7  <  1 


_^,Y/,.  xa/,.  \P   ...7/,,  \a/i,  \P/„.2  ,2,,2 
(12b)   T  = 


-a)^(k^)"(ky)P  ^  a3'y(k^)"(ky)^(co^-  ay) 
(kx'b'kx.'t)A  (kxb.k^b)D 
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(12c) 


T  =  T^  +  Tj.. 


All  the  transforms  in  (10)  with  denominator  AD  are  first 
derivatives  of  the  functions  whose  transforms  are  considered 
in  (12a).   We  consider  first  T^,  the  simpler  of  the  two  terms. 
For  7=1  we  may  write 


T  =  ^'-^K    (-i  l-)'^(-l  I-) 


P 


dCD 


J 


J 


(dk) 


^i(kT-cot ) 
'  2^2 


2(k^b^-(k.b)^)  /  CD  -  k^A^ 


"^  +  Vo 


a  +  3  =  1 


For  7  =  0,  so  that  a  >  1,  we  may  write 

r 


a-1 


P 


(dk) 


^i(k-r-cDt) 


2A~(??^(k^)/co 


Vo 


^  +  ^x^o 


a  +  P  =  2. 


Thus  it  suffices  to  consider 


S"  = 


:^  ct'  ^  p'    . 


r 


J 


i(k.r-a)t ) 

dcD|(dk)  -2 ^  ^  .  ^ 

((k.b)2-(I^.b^^) 


00  - 


^x^o 


^  -^  Vo 


a'  +  P'  =  1 
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We  note  that 


^'(|r)'-^(l+^')(ty)']  S-  =  H(t)(|^)°'(|^)^'|6(x-A^t)6(y) 


t   5(x+AQt)6(y)f  . 


An  elementary  calculation  gives 


+ 
S~  is  accordingly  a  rational  homogeneous  function  of  degree 

-1  in  x,y,t.  We  need  not  evaluate  S"  explicitly  as  a 

general  argument  will  show  it  does  not  appear  in  the  answer. 

On  the  other  hand  Tj.   satisfies  the  following  equation  : 


U  |f)'(-i  |j)°(-i  |y)'((i  If)' 


a  +  P  +  Y  =  2 
a  +  7  >  1 


Just  as  in  Section  III  of  I,  if  we  determine  Tj^  from  its 
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Fourier  transform  and  then  do  the  \k\      and  03  integrals  we 
obtain 

(14)   T,  =  ^  Ue    ;  ^  y ^ ^ 2 . 

Again  we  come  to  the  problem  of  the  interpretation  of  singular 
integrals.   The  finite  part  interpretation  of  the  integral  is 
appealing,  but  if  we  were  to  go  through  the  plane  wave  techniques 
suggested  in  I,  we  would  learn  that  again  special  irregularities 
associated  with  the  multiple  roots  of  D  might  appear.   But 
for  this  system  (characterized  by  a-1-7  >  1)  the  irregularities 
do  not  occur  and  we  may  write,  with  the  notation  of  I: 


(15)  T,  =  ^ 


L  -  ~^ 


At  the  points  €^(k^)^+  (l+€^)(k  )^  =  0,   or  at 

X  y 

2       2 
cos  6   =(l+e  ),  the  second  factor  in  the  denominator  vanishes. 

Let  C  be  a  curve  oriented  in  a  clockwise  sense  that  surrounds 

2      2 
the  values  cos  0  =  1+a   in  the  upper  half  plane  in  the  strip 

0  <  Re  e  <  TT     and  let  C  surround  no  root  of  D,  unless  it 
coincides  with  the  selected  values.   Let   C*'   be  a  similar  con- 
tour in  the  lower  half  plane  but  oriented  in  a  counter  clockwise 
sense.  Thus  we  have 
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r 


(16)  T^  =  aM 

kv   t 


d0 


r 


a 


-  ^ 


^k.rj,,k.?.2   _2, 


de 


(k,)  (k^)  (^)  ((^)  -  aj) 


rk-r 


2N,r  x2. 


C+C 


,  S(^,  k)(c^(ft j-^-Ki+c^Xk,,)^) 


X' 


y 


On  performing  the  second  Integration  we  obtain  a  rational 
homogeneous  function  of  degree  -1  in  x,y,t.   The  first  term, 
as  we  shall  see  belo'.\,  vanishes  outside  some  bounded  region 
of  the  x/t ,   y/t  plane.   We  expect  T  =  T.+  T^  =  Tj+   C  S~ , 
with  C~  some  constants,  the  fundamental  solution  of  a  hyper- 
bolic partial  differential  equation,  to  vanish  identically 
outside  some  bounded  set.   Thus  T»   and  the  second  integral 
in  (16)  must  cancel  in  this  region,  but  then  they  must  cancel 
everywhere  so  that  we  may  write  finally: 


(17)  T  =  ^W- 
4t  t 


a 


^Tr.^y     ^^-  ^  2 


(c.Oh-(c..c.)  d(^,  k^)(.2(^,)^(i..^)(kp2) 


We  might  readily  verify  that  (17)  is  just  the  result  obtained 
from  a  finite  part  interpretation  of  the  formal  Fourier  inver- 
sion of  (12a)  by  the  methods  of  Section  III  of  I,  but  we  shall 
leave  the  answer  in  its  present  form. 

In  I  the  change  of  variable  z  =  e^^^^"-^^^  transformed  the 
integral  1(12)  into  an  integral  of  a  rational  function  around 
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a  c 


losed  contour,  1(13), (1^),  expressible  in  terms  of  the 


residues  at  the  poles  enclosed  by  the  contour  (the  poles  occuring 
only  at  the  complex  roots  of  D(^^,  k))   1(22),   Here  too  the 
transformation  z  =  e^^^^~^^      turns  (1?)  into  the  integral 


around  a  closed  contour  of  a  rational  function  whose  only  poles 
inside  the  contour  occur  at  the  complex  roots  of  D(— ^,  k). 
We  might  write  explicit  forms  for  the  integral  analogous  to 
1(22),  but  as  they  are  not  particularly  illuminating  we  shall 
omit  them.   We  define  the  value  of  (1?)  to  be  H^j^(x,y,t), 
a  =  2,  p  =  Y  =  0;   H^y(x,y,t),   a  =  P  =  1,  Y  =  0|   Hyy(x,y,t), 
P  =  2,  a  =  Y  =  0. 

For  (17)  we  may  carry  over  essentially  all  the  general 
discussion  of  Section  IV  of  I  always  with  the  modification 
implied  by  (11)  included.   Thus  we  have  as  at  least  part  of 
the  ray  cone,  the  surface  given  by  1(19),  (20)  (with  the 
substitution  (11)).   Again  (see  Figure  1)  a  section  of  the 
ray  cone  consists  of  an  outer  region  A  which  the  disturbance 
has  not  yet  reached,  the  region  B  of  the  disturbance,  and  the 
lacunas  C.  There  is,  however,  an  addition  to  the  ray  cone 
coming  from  the  singularity  produced  as  the  roots  of 
D(^^,  ic)  approach  those  of  €^k^  +  (l+€^)k^.   A  direct 
calculation  (depending  on  the  fact  that   ( (^^^^)  -  a^)  occurs 

in  the  numerator  in  (1?))  shows  that  the  only  relevant  con- 

P    0  0 
fluence  of  roots  occurs  at  x  =  At,  y  =  0,   the  position 

of  a  pure  Alfven  wave.  Thus  there  are  added  singularity  points 

in  Figure  1,  the  points  E,E'   where  a  pure  Alfven  wave  would 
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be.   The  Integral  (I7)  has  the  character  of  a  compound  of  an 
Alfven  and  a  longitudinal  wave.  Just  as  It  should.   A  more 
complete  analysis  would  also  show  that  the  singularity  in  a 
derivative  of  T  near  E  or  E'  is  the  same  as  a  pure  Alfven 
wave,  i.e.,  it  is  of  the  form  6(x  ±  At)6(y). 

We  now  give  the  explicit  forms  for  the  physical  parts  of 
the  Green's  function: 

(18a)  p  =  (-  |^)Spj{{|)'-  a2)g^  +  f?  V  S  -yy 

V 

-   Ix  ^-^o  Su  3xx  ^  (V  X  %)-fl   A^cCf  G^+  f  G^p 
*   "*o(|  Oxy  +  I  Oyy) 

(i8i>)  (t.t)  =  (|j  Sp-  \^i$.t^m(^  -  4)3^*  sp  5;^+  2^  Oyy) 

.  a2(^.^J.[(|  V  F  «yy)"  ^  ^^(|  ^V  I  «xy) 
+  ao*o(^>^^'-(-^"°xx-  "  °xy) 
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(18c)     (B-b)  =  -(|f(t„-bHC^-^)(+  tF"[((f''-  ^o)°xx 

.  2p  0,,  .  (|)\^j  .   (S,.7){(b^(|)'-  a„2)G,, 

-  ^  3xy  -  ^  a„j-  S^  a^(^.#)(|  G,^.  f  5;j 

-  3pa,^jc^  1,(1  ^^  .  ^  S'„)-(i-^)|y(|  V  I  =yy> 


(184)      (u-T)  =  i.Al{t.^<%)-  |7(Su-J))(((|)^-  a2)G 


XX 


-fF^^xy-f^y-^(3u>-|r3p)(|aV^V) 


(I8e)      (^.^)   =     H(t)   ^  [((S^-n)-  A^Sg  )5(x-A^t) 


-^  t  C^xy)  +  ^^(^-V-S3)H^^  +eb^(n.7XS3)H^y 


-   ^(^>^^-bX^J(f  H^+f  H^y) 
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-fxci      ,    f-^    -is'b    \l,fX',^     „2x:r        .    2xy 


(I8f)      (B.J)    =   M%-i%  .   (Su-t)fj)    ((f)    -   a„-)G,,  +  ^  G,y 


^(f)\y)-^o(fFSp-S„     §j)5,^, 


'o  1^   (Sg-ixf  G,^.  f  Gyy).(^.^xSg)(|  H,,+  I  H^y). 


+  e(b-7^  S    )H 


(I8g)      (B-n)    =  H(t)    5M     (A^Sg   -   S^.n)5(x-A^t^. 

o   /  z 


Again  all  derivatives  operate  on  all  fvmctlons  to  the  rij^ht. 

With  reference  to  Figure  1,  all  the  functions  G,   H  vanish 
in  regions  A  and  C  and  become  infinite  as  one  approaches 
the  boundary  of  B.   On  the  boundary  of  B  the  singularity 
is  of  highest  order  at  D  and  D'.   For  the  functions  G 
there  are  no  other  singularities.   The  functions  H  are  also 
singular  at  E,  E',   the  position  of  the  Alfven  wave,  and 
behave  there  like  a  pure  Alfven  wave. 
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FOOTNOTES 

1.  H.  Weltzner,  "On  the  Green's  Function  for  Two  Dimensional 
Magnetohydrodynamlc  Waves  I",  NYO-2886,  New  York  University 
April,  i960.   We  shall  refer  to  this  report  as  I  and  use 
the  same  notation  here. 

2.  A  preliminary  account  of  these  results  was  presented  at  the 
November  I90O  meeting  of  the  Division  of  Fluid  Dynamics  of 
the  American  Physical  Society  at  Baltimore,  Maryland. 

1) .      See  A.  A.  Blank"  and  H.  Grad, "Notes  on  Magnetohydrodynamics 

VII  -  Fluid  Dynamic  Analogies",  NYO-6486,  Nev/  York  University 
July,  1958. 

'I.   H.  Grad  in  The  Magnetodynamics  of  Conducting  Fluids,  ed., 
D.  Bershader,  Stanford  University  Press,  1959- 
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OF=OF'=A 
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oH-y^, 
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FIG.  I 

SECTION    OF     THE   RAY    CONE    FOR  A  TWO 
DIMENSIONAL    MAGNETOHYDRODYNAMIC    WAVE 
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